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AN APPLICATION OF GROUP THEORY TO HYDRODYNAMICS* 

BY 

E. J. WILCZYNSKI 

It has been observed by Sophtjs Lie that the stationary motion of a fluid 
can serve as a perfect picture of a one-parameter group in three variables. So 
far as I know, neither he nor any of his followers utilized this fact for the pur- 
poses of hydrodynamics. It is the purpose of the present paper to do this. One 
of the advantages gained for hydrodynamics by this standpoint lies in the general 
conception. But another advantage is, as is always the case when a class of 
problems is investigated from a new standpoint, that from the group-theoreticaJ 
point of view, certain special cases are of exceptional interest, simplicity, and im 
portance, cases which otherwise would appear difficult and unpromising. 

§1. Relation between the steady motion of a fluid aud the theory 
of one-parameter groups. 

From the definition of the steady flow of a fluid, and Lie's conception of a one- 
parameter group, the following theorems will be seen to be true. 

1. Let X, y, zbe the Cartesian coordinates of any point of the fluid at the 
time t , and let a, b , c be the coordinates of the same material point at the 
time t = . If 

(1) x = ^{a, b, c; t), y = '>^{a, b, c; t), z = x{a, b , c; t) 

are the equations which describe the steady motion of any point of the fluid, 
they are also the flnite equations of a one-parameter group, the parameter 
being t . 

2. Conversely, if equations (1) are the equations of a one-parameter group, 
they can also be interpreted as describing the steady flow of a fluid. 

We shall be mainly concerned with the second theorem. Let 

be the infinitesimal transformation of a one-parameter group G , where u, v , 
w are arbitrary functions of x, y , z. The finite equations of the group are 
found by integi'ating 

* Presented to the Society (Chicago) December 29, 1899, under a slightly difCerent title. Ee- 
ceived for publication January 13, 1900. 
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(3) ^ = M(a5, y, »), ~^ = v{x,y,z), ~^^=w{x, y , z), 

with the condition x= a, y = h, z= c ior t=0 , if a, h , c are the variables 
which are transformed into x, y , a by the operations of the group. Thus u , 
V , w are the values of the velocity-components at the point x , y , z . 

Let X, Y , Z denote the components of the force acting upon the unit of 
mass, p the pressure, and p the density at the point x,y ,z. Then Euler's 
equations assume the form : 

(4) K. = X-^t-, K«=l--1^, K«,= Z-1^^, 
^ ' p ox p dy p dz 

and the equation of continuity becomes : 

(4«) ^P + P\Yx + Fy+-8.)-'- 

Moreover^ and p are connected by an equation : 

(46) P=f{p)^ 

where the nature of the function /"(p) depends upon the character of the fluid 
considered. 

Suppose that u , v , w are arbitrarily assigned functions of x , y , z . Then 
by integrating (3) we get x , y , z as functions of a , & , c , and t. The equations 
(1) , thus found, represent any one-parameter group G in space. We can find 
the physical conditions corresponding to every group G . First p may be found 
by integrating (4a). The simplest way to do this is to consider p as a function 
of a, b , c and t rather than of x, y, z. For then, since the motion is steady, 
i. e., since dp/dt = , Kp is simply dp/dt , and therefore (4a) may be written ; 

1 dp ( du dv dw\ 

whence 

where p^ denotes the value oi piov t = . In the case of steady motion, p^ is 
not an arbitrary function of a , 6 , c over any region of space. For if it were, 
p would be the same function of a; , y , z that /j^ is of a , 6 , c and therefore p 
would be an arbitrary function oi x, y, z which according to (4a) is obviously 
not the case. The reason for this difficulty is that although we have assumed 
explicitly that u, v , w are functions of x, y , z alone, we have as yet made no 
explicit assumption of the same character about p , equation (5) being universally 
true whether the motion be stationary or not. 
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If /S be a surface, cutting all of the lines of flow once and only once, p can 
obviously be assigned arbitrary values for all points upon this surface. Of 
course S may consist of several separate pieces. When any element of the 
fluid reaches this surface, its density must then, if the motion is stationary, as- 
sume the value corresponding to that point of S. Thus and only thus can the 
motion be stationary. 

After p is obtained (46) gives p , and from (4), ^, Y, Z can be found. If 
u, V , w are real, and p^ is positive, p from (5) is positive as it should be. If 
p should become negative, the fluid would tear. 

In the case of an incompressible fluid, these results would be somewhat dif- 
ferent : p is not then a function of p . We then have 

du dv dw 

so that u, V, w must be taken subject to this condition. The surfaces 
p = constant are made up of lines of flow, for the density of an element of an 
incompressible fluid can not change. From one line of flow to another the density 
can vary according to an arbitrary law, for (4a) has reduced to Kp = , which 
says only that the surfaces p = constant are invariant under all transformations 
of the group, i. e., that they are made up of lines of flow, which we have already 
assumed to be the case. Since p is not in general a function oi p, in this case 
p can be assumed to be an arbitrary function of x, y , z and then ^, Y, Z 
are determined from (4). 

If the fluid is homogeneous, without being incompressible, p = constant, and 
then from (4a), since Kp = , we get 

■du dv dw 

i. e., although the fluid may not be incompressible, if it is homogeneous and has 
a steady motion, it moves just as an incompressible fluid of like density would. 
But if the fluid is not incompressible, ^ is a certain function of p , and since 
p = constant,^ = constant all through the fluid. The forces X, T, Z are 
again obtained from (4). 

The one-parameter group, generated by K/", has two families of invariant 
surfaces, viz., the cxs ^ surfaces : 

tij(x , y, z) = constant, u^^x , y , z)= constant, 

where ti^ , u^ are two independent solutions of Kf= . The intersections of 
these surfaces are the lines of flow of the fluid. 

Like all one-parameter groups, ours can be written in the form : 

(7) n.{x,y,z) = n.{a,b,c){i = l,2), W{x, y, z) = W{a ,h , c) + t . 
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Thus all steady fluid motions can he put into this form, and every set of 
equations of this hind represents a steady fluid motion. 

It should be remarked, however, that it is not necessary that the functions li^ 
and W should be represented by the same formula for all portions of space. 
Thus, we have more generally 

(8) 

TF,(a3,y,z)= Tf,(a,6,c) + < (&=l,2, 3, •• •), 

the formulae with index h being valid only for a certain region J?^ of space. 

Thus while the point a; , y , a is in R^, fl,^ and W.^ are used, when it is in 
M^ , the next adjoining region, il^^ and W^ are employed, etc. Of course, for 
continuity, the values oi x, y , z obtained from the first set of equations must 
coincide with those found from the second set for all points upon the boundary 
between J?j and R^, \. e.,x,y ,z must be continuous functions of a , 6 , c , t . 
The same is not necessary of m , v , to , i. e., in general the direction of the 
motion will undergo a discontinuous change on the bounding surface between 
two such regions. Such discontinuities will actually occur, for instance, if we 
consider the flow of a fluid in a channel whose walls are portions of analytical 
surfaces, which meet at angles different from 180°. 

These considerations go a great way towards the solution of the general 
problem of the steady flow of a fluid in a channel bounded by any arbitrary sur- 
faces whatsoever. An infinity of equations of the form (8) can be written 
down in every case, which fulfil the boundary conditions, and which represent a 
possible fluid motion. Among all these, it will be necessary to pick out those 
which correspond to the external forces at work, i. e., which give the values to 
Ji. , Y , Z belonging to that particular problem. 

In most cases the forces X, Y, Z have a potential, so that 

^^) ^=-^' ^^'dy^ ■^=fe- 

In all such cases u, v , w cannot be chosen at will, but are subject to the 
condition, that the values of ^, Y, Z computed from (4) can be written in 
the form (9). If we put 

(10) P = /* 

equations (4) become : 

^1^ v S{V-P) ^ d(V-F) ^ 8iV-P) 
and therefore, if the forces X, Y, Z have a potential, 
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(12) Kudx + Kvdy + Kwdz = d{V-P) 

must be a complete differential. This is the condition which u, v, w must 
verify in this case, or, what amounts to the same thing, 

dKv dKw ^Kw dKu _ dKu 8Kv _ 

^ ' dz dy ~ '' dx dz ~ '' dy dx ~ ' 

§ 2. The fluid motion expressed by the general projective group. 

We shall confine ourselves to the case in which the forces have a potential 
and we shall find that the most general ternary projective group, which can ex- 
press the Steady motion of a fluid, if the forces have a potential, is a linear 
group. 

For convenience in this paragraph we shall write x^, x^, x ^, Mj , u^, u^ in- 
stead of X , y , z , u , V , w. 

The most general projective infinitesimal transformation is a linear combina- 
tion of the following fifteen : 

^•=£' "'A=^a' ^ii.^,Pu=^i (f,fc = l,2,3). 
We can write it : 

(14) Kf=±u,^, 
where 

3 X 

(15) Mi = 7.- + E <:A + a'X <^,^K (i = 1 , 2 , 3). 

k=l k=l 

We must find the necessary and sufficient conditions in order that 

3 

5] Kuxix. 



may be a complete differential. 
We obtain 

3 

(16) v^ = Ku. = (cjX, + c^x, + CjXj)!*. + 23 (c^ + c,a;>, (i = 1 , 2 , 3). 

But if the forces have a potential, equations (13) of § 1 must be verified, the 
first of which is 



We find 



^3 _ ^2 _ Q 

dx„ dx. 



i:((c..3+c,x3)||;-(e..+c,x,)g;) 



dv, dv. 

3 
+ 
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Now, from (2), we find 
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du 

i 

dx. 



- = c^i + c/c. + B^^c,x^, 



3 

c 

K=l 



where 8^. = for i'^j , and S^^ = 1. 
Thus we obtain : 

^ - ^ = c,«3 - C3«2 + [^23 - C32 + (c^aJj - CjXj)] gc;,a;;, 

3 3 3 

*=1 A=l A=l 

When this is arranged according to powers of x^, x^, x^, the coefficient of 
each term must vanish. This gives the following equations, obtained by put- 
ting the coefficients of the quadratic terms equal to zero : 

which gives at once c^ = Cj = . Using one of the other equations (13) §1, we 
should find Cj = also. 

Equating the other coefficients in {dvjdx^ — {pvjdx^ to zero, we get only one 
other condition, and thus we have proved the theorem : 

An arbitrary projective one-parameter group cannot represent the steady 
motion of a fluid under the influcTice of forces possessing a potential. To 
represent such a motion, the group must he a linear group. If 



(17) 

where 

(18) 






h = % + Z hi^k ' 



3 

*=1 



is the infinitesimal transformation of this linear group, the coefficients c^^. must 
verify the three conditions : 



(19) 



3 

3 
Z (<'*2<'l* - Cj,C2,) = , 

which necessary conditions are also sufficient. 
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§ 3. Fluid motion expressed by the ternary linear group. 

In equations (18) we can put 7^ = without loss of generality. For the 
motions in which 7; 4= ^ only differ from those for which 7^ = by a translation. 

We assume then that the origin remains fixed during the entire motion, and 
we will resume the less symmetrical but more convenient notation, x, y , zi u, 
V , w instead oix^, x^, x^; u^, u^, u^. Let 

■ u = ajX + b^y + CjZ; 

■ v = a^ + h^+c^, 

,w = afi + b^ + c^. 

The components of the velocity of rotation of any element of the fluid are : 



(20) 



(du dw\ 



so that f , »7 , ? are the same for all points of the fluid. Hence at all points 
the axis of rotation has the same direction, and the angular velocity of rotation 
is the same, viz.: 

Let us take our z axis parallel to this direction. Then 

1=^ = 0, ?=w, 
and hence, for this choice of coordinates, 

*3 = <'2 ' Cj = 083 , a^ = 6, + 2ft) , 
so that we can write more simply 

M = Ax + {H — (o^y + Gz, 

(21) ■ v={H+m)x+ By + Fz, 

w= Gx+ Fy+ Cz, 

If the forces have a potential, equations (19) give 

{A + B)m=Q, ii^ft) = , (?« = . 
Therefore, either o) = , or 

^ + 5 = 0, i^=0, G = Q. 
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If ft) = , we have a velocity potential 

<j> = ^ {Ay? + By" + <7z^ + 2Fyz + Wm + 2Hxy), 

so that the lines of flow are the orthogonal trajectories of the system of similar 
surfaces of the second order : 

4> = constant. 

In the second case (21) reduces to 

u = Ax + (ff— <>>)y , 

(22) ■ v = {H-\-m)x— Ay, 

w = Cz . 

The motion can be decomposed into a rotation and a motion with a velocity 
potential. But it is simpler to treat it without so decomposing it. This case 
again separates into two distinct sub-cases. 
In the first sub-case, 

A'+II'-CO:^0. 

Then, if 



p^= + ^A"+H"-m = p, p^=-p^=-p, 
H- m A + p. 






(23) ^ ^ _ 

A-p^ H+m 

the finite equations of the group can be written as follows : 

(24) x — X^y= e-i'\a — Xfi) , x—\^ = e+i''(a — Xjb) , z—ce 

where a, b, c are the values of x, y, z for t = , 
The surfaces : 

z-i'(x — \y)-^= constant, zf^x — \^)~^= constant, 

which are real or imaginary cylinders, are therefore invariant. Their real in- 
tersections are the lines of flow. 
The function 

(x - \y) {x - \0) 

is also invariant. Therefore the projections of the lines of flow upon the xy 
plane are the similar conies : 

.„., , 2A w — H . 

(25) XT — -rr—. — '- as/ -I XT- V = coustaut , 

i. e., the lines of flow are situated upon a cylinder of the second order, whose 
elements are parallel to the z axis. 
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The physical conditions, corresponding to this motion, are obtained from the 
considerations of § 1. We have, in particular, from (11), 

(26) F- p = J 1{A^ + H^- m){x^ + f)-¥ cV] + constant. 

The motion, here considered, becomes periodic if p is purely imaginary, and 
c = , as may be seen by solving (24) for x , y , z . The lines of flow are 
then similar ellipses in parallel planes. It can be shown that this motion is 
possible for an incompressible homogeneous fluid filling an ellipsoidal space, the 
particles attracting each other according to Newton's law. This has been 
studied in detail by Dirichlet and Dedekind.* It can also be easily shown 
that the analogous result is true for an elliptic cylinder of infinite length, pro- 
vided that H and a satisfy the equation : 

where p is the density, andy a numerical constant, depending upon the units 
employed. If a and h' are the semi-axes of the elliptic section of the cylinder 
their ratio will be obtained from the equation : 






-H 

In the second sub-case, 



H-(o A 



A H+ ft) ■ 
The finite equations of the group can be written : 

If we chose the plane a; — Xy = as ass-plane, and the plane Xx + y = as yz- 
plane, the planes y = constant are invariant, and the lines of flow are the curves : 

z = constant e^" * , y = constant , 
i. e., exponential curves in these planes. 

§ 4. Fluid motion expressed hy a linearoid group. 
Let us consider next a one-parameter group, generated by 



*Crelle's Journal, vol. 58. 
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•/•i » ■ ■ ■ ' Xs heing functions of z only. It will be of the class called linearoid by 
me in another connection. Then, 



Kw = (x,<^, + X2^i + x'iXy + (XA + Xifi + X3X2)y + Xi<^3 + Xaf 3 

+ X3X3 + xixi*' + (xix2 + X2X>y + X2X22/2 » 



(28) 



where the accents denote derivatives with respect to z . 
Write this as follows : 



(29) 



K« == \x^ 4. 2fi^xy + vy + 2p^x + 2a-^y + 2tj , 
Kv = X^2 ^ 2/t^y + vy + 2p^ + 20-2^ + 2t, , 
Kw = XjCC^ + 2/t3a;y + v^ + 2/0303 + 2<7-3y + 2X3 ; 

the conditions for the integrability of 

Ku dx + 'Kvdy + Kto • dz , 

i. e., for the existence of a potential for the corresponding fluid motion, become ; 

>, - Xj = , I/, — ^2 = , o-j - P2 = , 
P2 - M3 = , 0-; - 1/3 = , t; - 0-3 = , 

p[-\i=Q, c.;-,t.3=o, t;-,,3 = o, 

^; = /^i = "1 = '^a = 1^2 = ''2 = . 

2(^1 - \) = (<^; - "-nxx + <i>:x, = , 
sc"! - ^2) = - ^; X. + (2</>; - f[) X2 = , 



(30) 

Hence 

(31) 



so that we have two cases to distinguish, according as Xi = Xi — ^ ^ ^^ ^^^ deter- 
minant 

(32) (,^; - 2^;)(2</.; _ ^;) + ,^;^; = o . 
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Let US consider first case I, in which ^^^ '^^= Q . Then we find from the 
definition of these quantities, 

\^ = /ij = l/j = \^ = ^^ = V^ = \3 = ^3 = 1.3 = P3 = 0-3 = , 

and from (30), 

Pi = P2 = ""I = '^2 = "^i = < = *^ ' 

SO that pj , p^, o-j , o-j , T, , Tj are constants. 
Now since ^^ = ')^^= , we have 



(33) 



(2p, = <t>l + <l>,f,, 2p, = ^,(<^,+ f,), 
I 2cr. = U<f>i + ^2) . 2cr, = <l>,f^ + fl . 
The only further condition is a-^ = p^, whence either 

^1 + f 2 = , or 4>2='^i- 

"We will first assume that <^j + i/r^ = , and call this case la. Then we 
have, in this case, 

<^. = P2 = 0, <^J + <^,^i=2p, = 2<r,= 2.r, 

where o- is a constant. Thus in case la the functions Xi ^nd X2 vanish, and 
between the other functions we have the relations : 



(34) 



<t>i +^^2= ^' ^l + ^2"*^! == 2a- == constant, 

't>A + <^2t3 + ^S'i>'s = 2Ti , 

fA + tifs + Xs^s = 2t2 , 



where r^ , t^ are also constants, and where <f>^, <f>.^, x^ may be taken as arbi- 
trary functions of z . 

In the second case, case 15, we have 

^2 = "f 1 ' 
whence 

where (/>, and sjr^ may be taken to be arbitrary functions of z , and the equa- 
tions for ^3 , i/tj , Xs are of the same form as in (34). 

If the fluid is incompressible we have further, X3 = constant in case la, and 
in case 16, we find 



Xi= — j i4>i + -^^dz + constant, 
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as is seen at once from the condition of incompressibility, viz. : 

du dv dw 

dx dy dz~ ' 

"We shall not at present farther discuss the fluid motions thus analytically de- 
termined. The group theoretical considerations have enabled us to discover 
a considerable class of fluid motions, whose entire theory depends only on linear 
differential equations. 

Let us now carry out the corresponding investigation for case II, in which Xi 
and X2 do not both vanish. We have then 

(32) (,^;_2Vr;)(2^;-^;) + ,^;t; = o, 

and further, as before (equations (30)) Xj , 11^, v^, \, fi.^, v^ are constants. 
By the definition of these quantities, we have 

\f[Xl=\^ f2Xl + V^iA:2 = 2M2, f2Xi==''2' 

Therefore 

(37) "^f + ^^^-a^XX^^O (»=1,2), 
and hence 

(38) \:2(l,:v, = X,^.2^i,:v,= a^.0^.1, 

where a and denote the ratios \j : Vj and 2/*, : y, . 

Multiplying (32) by xlxl » ^^^ using (36) and (38), we obtain a third quad- 
ratic for Xi and X2 » ^2. : 

(39) (^iXi-'^^^/^O, 
so that also 

XI : 2X^v^ : vl = a: /8 : 1. 

Now if we rotate the x axis through an angle 6 , the form of w , u , w is not 
changed, and w becomes 

Xi{x cosd—y sin 0) + X2 {^ sin + y cos 0) + Xs 

= (Xi cos ^ -f X2 sin (?) K -f- ( - Xi sin ^ + X2 ^'os 0)y + Xi- 
If we put 

V, 

tan ^ = ^ , 
according to (39), the coefficient of y in the transformed w will vanish. 
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We can therefore assume from the beginning, without loss of generality, that 
X^ = . Having already treated the ease x^ = j^^^ , we will now assume 

From (36) we have, in this case, 

Al Al 

Since Vj = , (30) gives /tj = . Therefore, if 






we have, since /tj = X^ , 

J <^j = \<^ + a , <^, = 2\4> + fi , 

where a , y8 , 7 , h denote four arbitrary constants. 
We find further , 

\ = X'xXi » Ms = ' "3=0- 

Therefore, from (30), we notice that p^, o-^, a^ are constants. Now we have, 
with the values for </)j , ^^ , i/tj , -<|t-^ just found, 

2^, = 2\^</.^ + \(yS + 27).^ + /37 + 8% 

so that X2 nmst vanish. Then 

2^, = y8(V + a + S), 

so that \ must also vanish. From the definition of p^ , which must be a con- 
stant, follows : 

(41N 2p,-7(a + g) _ (^-7)(a+a) 

^ ' ' Xi Xi ' 

since p^= a-^. Similarly we find, from joj = X^ = x'lXi > 
,.^^ ,, e — a^ — (87 

(42) «^; = xi + — ^— ~^, 

Al 

where e denotes another constant. 

Substituting all these values into the only two equations (30) not yet used, 
viz. : 72 = 0-3, T, = P3 , we find the following conditions for Xi and ^3 : 
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,y(g _ g^ _ ^i) + h{^ -i){a + g) 



(43) 



(7 - ^)X. + -, 



+ 



(^-^x« + ^)i(|)=«' 



a(e-a^-/g^) + ^(^-^)(a-g) f 



+ (e-a^-^,)|(|)-X,x;=0. 



The fluid motion, in this case, is therefore given by the equations : 

■ M = ax + /8y + <^3 , 

(44) I -y = 7X + 5?/ + ■f 3 , 

. w = XiX + X3 , 

where a , /8 , 7 , S are constants, where Xi ^^^ X3 ^^^ obtained as functions of 
z by integrating (43) , and where ^3 and i/r^ are determined from (41) and (42). 
If the fluid is incompressible, we have 

du dv dw r. / ' n 

i. e., 

Xi = V, X3 = - (a + S)z + ^ , 

where rj and d are constants. Moreover d can be put equal to zero, if the ccy- 
plane be appropriately chosen. Substituting these values in (43), we get rela- 
tions between r] and the other constants. But as will be seen from (41) and 
(42) the group in this case becomes a ternary linear one, all of the coefficients 
being constants. 

Thus Case II of this paragraph gives nothing new for incompressible fluids. 

The detailed discussion of the fluid motions here determined, as well as those 
of the preceding case, will be left for a future occasion. 

University of California, Berkeley. 



